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Abstract: We present a determination of the strong coupling constant from a fit of 
QCD predictions for six event-shape variables, calculated at next-to-next-to-leading order 
(NNLO) and matched to resummation in the next-to-leading-logarithmic approximation 
(NLLA) . These event shapes have been measured in e''"e~ annihilations at LEP, where the 
data we use have been collected by the ALEPH detector at centre-of-mass energies between 
91 and 206 GeV. Compared to purely fixed order NNLO fits, we observe that the central 
fit values are hardly affected, but the systematic uncertainty is larger because the NLLA 
part re-introduces relatively large uncertainties from scale variations. By combining the 
results for six event-shape variables and eight centre-of-mass energies, we find 

a^(Mz) = 0.1224 ± 0.0009 (stat) ± 0.0009 (exp) ± 0.0012 (had) ± 0.0035 (theo), 

which improves previously published measurements at NLO+NLLA. We also carry out 
a detailed investigation of hadronisation corrections, using a large set of Monte Carlo 
generator predictions. 

Keywords: QCD, Jets, LEP Physics, NLO and NNLO Computations, resummation, 
strong coupling constant. 
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1. Introduction 

Event-shape distributions in e"'"e'~ annihilation have been measured with high accuracy 
at LEP at centre-of-mass energies between 91 and 206 GeV [1-5] and at the SLAC-SLD 
experiment at 91 GeV [6], as well as at the DESY PETRA collider at lower energies, e.g. by 
the JADE experiment [7]. Event-shape obscrvablcs are infrared-safe variables designed to 
describe the structure of the hadronic final state. At leading order in perturbation theory, 
e+e" annihilation to hadrons occurs via e+e" qq and subsequent hadronisation to stable 
hadrons, resulting in a back-to-back (two-jet-like) structure of the event. At higher orders, 
gluon radiation off quarks will lead to deviations from this two-jet structure. 

Fixed-order QCD corrections to event-shape distributions were calculated some time 
ago at next-to-leading order (NLO) [8-10], and more recently at next-to-next-to-leading 
order (NNLO) [11-16] for the six event-shape observables thrust T [17] (respectively r = 
1 — T), heavy jet mass Mh [18], wide and total jet broadening By/ and Bt [19], C- 
parameter [20] and the two-to-three-jet transition parameter in the Durham algorithm, 
— Inya [21]. The definitions of these variables, which we denote collectively as y in the 
following, are summarised e.g. in [13,22]. 

As the fixed order expansion is reliable only if the event-shape variable is sufficiently 
far away from its two-jet limit, i.e. away from y — > 0. This is because large logarithmic 
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corrections spoil the convergence of the perturbative expansion in the two-jet region, in- 
dicating the sensitivity to multiple soft gluon radiation. To obtain a reliable theoretical 
prediction in the full kinematical range, it is therefore necessary to resum these logarithms 
to all orders in perturbation theory and to match the resummed result to the fixed order 
calculation. 

Until very recently, the theoretical state-of-the-art description of event-shape distri- 
butions over the full kinematic range was based on the matching of the next-to-leading- 
logarithmic approximation (NLLA) [23] onto the fixed next-to-leading order [8-10] calcu- 
lation. Now that the NNLO results are available, the matching of the resummed result 
in the next-to-leading-logarithmic approximation onto the NNLO calculation has been 
performed [24] for all six event shapes mentioned above in the so-called In i?-matching 
scheme [23]. It was found that the difference between NLLA+NNLO and NNLO is largely 
restricted to the two-jet region, while NLLA-I-NLO differed from pure NLO in normal- 
isation throughout the full kinematical range. For the thrust distribution, logarithmic 
corrections reaching beyond the NLL approximation have been calculated recently [25] 
using Soft-Collinear Effective Theory (SCET) [26]. 

In the theoretical description of event-shape observables within perturbative QCD, the 
only free parameter is the strong coupling constant a^, such that a fit of QCD predictions 
to the data for these observables lends itself to determine the strong coupling constant with 
high precision. Several determinations of based on NNLO results have been performed 
recently: using NNLO predictions [13] for the six event shapes listed above, a determination 
of ag based on ALEPH data [1] has been performed [27] , where the systematic uncertainty 
from renormalisation scale variations was found to be reduced by a factor of two as com- 
pared to the fit based on NLO predictions only. After the NNLO-I-NLLA calculations for 
the six event shapes have become available [24], a determination of based on JADE data 
has been carried out in [28]. 

Further, there are several studies based on thrust only: Using ALEPH and OPAL data 
for the thrust distribution and combining the theoretical NNLO prediction with infrared 
logarithms resummed within the SCET formalism, a precise determination of a^, has also 
been performed in [25]. A re-evaluation of the non-perturbative contribution to the thrust 
distribution is presented in [29], where the NNLO results of [11, 13] have been matched to 
resummation at NLL accuracy and then used for a combined determination of both the 
low-scale effective coupling uq and as(Mz). 

A very recent study of non-perturbative corrections based on moments of event shapes 
has been carried out in [30,31], ref. [31] containing also a determination of both cto and 
as{Mz)- However, these studies are based on NLO calculations only. NNLO predictions 
for event shape moments can be found in [32]. 

The agreement in the two-jet region of the matched prediction with hadron-level data 
is still far from being perfect, the discrepancy being attributed mainly to non-perturbative 
hadronisation corrections, but also to missing subleading logarithms, electroweak correc- 
tions and quark mass effects. In fact, based on the results of [25], one can estimate that 
the subleading logarithms can account for roughly half of the discrepancy between the 
parton-level matched NLLA-f-NNLO prediction and the hadron level data. 
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While significant progress has been made for the perturbative calculations, the non- 
perturbative corrections for hadronisation needed to extract the value of as from event- 
shape distributions are still obtained from Monte Carlo event generators based on leading- 
logarithmic (LL) parton showers and fragmentation models [1,27]. The hadronisation itself 
is presently parametrised by string- or cluster fragmentation models but the simulation of 
the multi-parton final state can now be performed at NLO+LL, which is in principle more 
consistent with the NNLO+NLLA calculation we use in our fits. We therefore investigate 
the performance of this type of generators taking HERWIG++ [33] as reference, which 
represents a modern event generator allowing optionally the inclusion of NLO calculations 
according to different schemes. 

In this paper, we present a determination of the strong coTipling constant based on 
matched NLLA+NNLO results for six event-shape variables. Hadronisation corrections 
and quark mass effects (at least to NLO [34]) have been included in the procedure. We first 
review the theoretical framework in section 2 before proceeding to the comparison with data 
in section 3. The method used for the fit follows closely that described in [1,27], but some 
improvements to the method used in [1] have been made which will be explained in section 
3. Systematic uncertainties for individual determinations of (Q) from different variables 
at different energies are presented in section 4. Combined results are given in section 5 
and further systematic studies concerning the hadronisation corrections, the scheme of 
matching NLLA to NNLO, the normalisation and quark mass correction procedure as well 
as the combination method are discussed in section 6. Finally our findings are summarized 
in section 7. 



2. Theoretical Framework 

The fixed-order QCD description of the experimentally measured event-shape distributions 

1 da 
o-had dy 

starts from the perturbative expansion 

— ^(y,Q,/i) = as{fi)^{y) + a'i{ii)^{y,x^) + al{fi)^{y,x^,) + 0{ai) , (2.1) 
(To ay ay ay ay 

where 

" 27r ' '^'^ " g ' 

and where A, B and C arc the perturbatively calculated coefficients [13] at LO, NLO 
and NNLO. They have been computed with the parton-level event generator program 
EERAD3, which contains the relevant matrix elements with up to five external partons [35- 
38], combined using an infrared antenna subtraction method [39]. A recently discovered 
inconsistency in the treatment of large-angle soft radiation [14] in the original EERAD3 
implementation has been corrected, resulting in numerically minor changes to the NNLO 
coefficient in the kinematical region relevant to the phenomenological studies here. In 
the deep two-jet region, e.g. (1 — T) <C 0.02, these soft correction terms turn out to be 
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numerically significant. They account for an initially observed discrepancy between the 
EERAD3 results and the logarithmic contributions (computed within SCET) to the thrust 
distribution to NNLO [25], which are now in full agreement. 

All coefficients are normalised to the tree-level cross section for e+e~ qq, ctq. For 
massless quarks, this normalisation cancels all electroweaJc coupling factors, and the depen- 
dence of (2.1) on the collision energy is only through Og and x^. Summation over massless 
quark flavours in o"o and da/dy results in a constant normalisation factor which cancels 
exactly in the ratio of these quantities. 

Predictions for the experimentally measured event-shape distributions are then ob- 
tained by normalising to (Thad as 

' ?(y,Q,/^) = ^?^-?(y,Q,M). (2.2) 



o-had dy ' ' o-had(Q,M) o-o dy 

For massless quarks, all electroweak coupling factors cancel in <To/c7had- 

In all expressions, the scale dependence of is determined according to the three-loop 
running: 

where L = 2 ln(/x/A2^^)) and A are the MS-scheme coefficients listed in [13]. 

The assumption of vanishing quark masses, which was used in all expressions for dif- 
ferential distributions up to here, is only partly justified, especially for the LEPl data, 
where bottom quark mass corrections can be relevant at the per cent level [40]. The effect 
scales with M^/Q^ and decreases to 0.2-0.3% at 200 GcV. We take into account bottom 
mass effects by retaining the massless Np = 5 expressions derived above and adding the 
difference between the massless and massive LO and NLO coefficients A and B [34], 

-^^{y,Q,t^) = ^ (ii-n{Q))^ +n{Q)^ ). (2.4) 

Chad Cly CThad \ O-V masslcss '^2/ massive/ 

A pole b-quark mass Mb = 4.5 GeV/c^ was used and Standard Model values were taken 
for the fraction rh{Q) of bb events. In this case, the electroweak coupling factors no longer 
cancel in the ratio ao/ahad, and the summation over quark flavours has to be carried out 
explicitly. 

For (Thad an NNLO calculation {0{al) in QCD) [41] including mass corrections for the 
b-quark up to 0{as), and including the leading mass terms to 0{a^), was used to calculate 
the correction cxo/cThad- A genuine 0(a5)-expression for (2.2) can be obtained by expanding 
the ratio ao/ahad, as done in [13]. 

Next-to-leading order electroweak corrections to event shape distributions in e+e" an- 
nihilation were computed very recently [42] . Using the event selection cuts and event shape 
definitions as applied in the experimental analysis [1], one observes substantial electroweak 
corrections to (Thad and da/dy. These corrections cancel to a very large extent in the ratio 
(2.2). In the kinematical range used in the determination below, the total electroweak 
corrections are at the level of two per cent at LEPl and at most five per cent at LEP2. 
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Genuine weak corrections from virtual massive gauge boson loops or fermion loops amount 
to one per mille or less at both LEPl and LEP2. The corrections are thus much lower than 
initially anticipated from partial calculations of higher-order electroweak contributions [43] . 
Since the experimental data were corrected for photon radiation effects using PYTHIA, it 
is not straightforward to include the electroweak corrections, and requires further study. 

The resummation of large logarithmic corrections in the y — ^ limit starts from the 
integrated cross section: 

RM,,)^J- f'Jfd^,., (2.5) 

0-had Jo dx 

which has the following fixed-order expansion: 

Riy,Q,,i) = 1+ A (y) a, (/x) + B (y, x^) o?, (/x) + C {y, x^) (m) + 0{ai) . (2.6) 

The fixed-order coefficients A, B, C can be obtained by integrating the distribution (2.1) 
normalised to cjhad (2-2) and using i?(ymax) Q, fJ-) = 1 to all orders, where ymax is the 
maximum kinematically allowed value for the event-shape variable y. 

In the limit y — one observes that the perturbative ^"-contribution to R{y) diverges 
like ctgL^", with L = —Iny {L = — ln(y/6) for y = C). This leading logarithmic (LL) 
behaviour is due to multiple soft gluon emission at higher orders, and the LL coefficients 
exponentiate, such that 

\n.R{y) ^ Lgi{asL) , 

where gi{asL) is a power series in its argument. 

For the event-shape observables considered here, and assuming massless quarks, leading 
and next-to-leading logarithmic (NLL) corrections can be resummed to all orders in the 
coupling constant, such that 

R {y, Q, Ai) = (1 + Cias) e(^9i(a.^^)+92(a.L)) ^ (2.7) 

where terms beyond NLL have been consistently omitted, and n = Q {x^ = 1) is used. By 
differentiating expression (2.7) with respect to y, one recovers the resummed differential 

event-shape distributions, which yield an accurate description for y ^ 0. 

Closed analytic forms for the LL and NLL resummation functions gi{asL), g2{asL) 
are available for r [44], Mh [45], Bw and Bt [46,47], C [48] and ys [49]. They can be 
expanded as a power series, such that 

oo i+1 

lni?(y, Q,IJ)=Y.11 Gi,,+2-n4^'+'~" • (2.8) 
1=1 n=l 

In order to obtain a reliable description of the event-shape distributions over a wide 
range in y, it is mandatory to combine fixed order and resummed predictions. However, in 
order to avoid the double counting of terms common to both, the two predictions have to be 
matched onto each other. A number of different matching procedures have been proposed 
in the literature and for a review we refer the reader to Ref. [22]. The most commonly 
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used procedure is the so-called In i?-matching [23]. In this particular scheme, all matching 
coefficients can be extracted analytically from the resummed calculation, while most other 
schemes require the numerical extraction of some of the matching coefficients from the 
distributions at fixed order. Since the fixed order calculations face numerical instabilities in 
the region y ^ 0, the numerical extraction of matching coefficients is prone to large errors. 
Therefore we restrict ourselves to the In i?-matching, studying two different variants of the 
latter for the present analysis. In the In i2-matching scheme, the NLLA+NNLO expression 
is 

In {R {y, as)) = L gi {asL) + 52 {ctsL) (2.9) 
+ as {A{y) - GiiL - G12L'') + 

+ a| {b (y) - (y) - G22L' - 

+ a| [c {y) -A{y)B (y) + ^-A^ {y) - GssL' - + 0{ai) . 

The matching coefficients appearing in this expression can be obtained from (2.8) and are 
given in [24], where we refer the reader to for more details. Quark mass corrections to this 
expression are included by retaining the mass dependence of the fixed-order coefficients 
A and B, which follow from the mass-corrected coefficient functions (2.4) and the mass- 
dependent total hadronic cross section. 

The coefficients in (2.9) explicitly depend on x^, thereby stabilising the scale depen- 
dence of the theoretical prediction: 

as ^ a,(/x) , (2.10) 



B{y) ^B{y,fi) = (3oynx^A{y)+B{y) , 

C{y)^C {y, li) = {(3o Inx^f A (y) + Inx^ [2 /3oB (y) +i3iA (y)] +C{y) , 

(2.11) 

92 (asL) 92 {asL,f/) = 52 {asL) + ^ {asLf g[ {agL) Inx^ , (2.12) 
G22 G22 (m) = G22 + /?oGi2lnx^ , 

G33 ^ G33 (Ai) = G33 + 2/3oG23lnx^ , (2.13) 

where g'l denotes the derivative of gi with respect to its argument. This scale variation also 
exemplifies a tension between NLLA and NNLO, since the NNLO coefficients compensate 
the renormalisation scale variation of ag up to two loops, while the NLLA coefficients only 
compensate the one-loop variation. A fully consistent matching, including the full scale de- 
pendence, is therefore only accomplished by combining NLLA-I-NLO or NNLLA-I-NNLO. 
In order to assess the effect of this incomplete compensation of scale-dependent terms, 
we have computed the two-loop terms proportional to in the above resummation and 
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matching functions, i.e. the scale-dependent logarithms appearing in and in the associ- 
ated matching coefficients G21 and G32, and recomputed the theoretical error in this new 
matching scheme, which we call the lni?(/Li)-scheme. In this scheme the NLLA-I-NNLO 
expression becomes 

In {R{y, as)) = L gi (a^L) + 52 ("s^) + "5 53 (asL) 
+ as {A{y)-GiiL-Gi2L^) + 

+ 4 (b iy) - _ c^iL - G22i^' - G23L^ 

+ 4 (^C (y) -Aiy)B (y) + ^A' (y) - G32L2 - ^33^^' - ^34^^") + O(a^) , 

(2.14) 

where the /it-dependence of 53 is 
93 {asL) 93 {asL, /x^) = c/3 (asL) + (asL) Inx^ Po92 (oisL) + ^ (asL) g[ (asL) 

4 («.L)^(l.x,f ^(A(L,a„sL))). (2.15) 

The "bare" 53 (asL) is not known and is put to zero, whereas the renormalisation scale 
dependence is proportional to the derivatives of gi and 52- In order to have a correct 
compensation of the renormalisation scale dependence the following further substitutions 
are made 

G21 G21 (/x^) = G21 + PoGulnx^, 

G32 ^ G32 = G32 + (/3olnx^)' G12 + Inx^ (2/3oG22 + 2/3iGi2) . (2.16) 

As discussed in detail in Section 6 below, the resulting theoretical error is consid- 
erably lower than in the standard In i?-matching scheme, and comparable to the error 
obtained in [25], where the thrust distribution was computed beyond NLLA and matched 
to NNLO [11,13]. 

In order to ensure the vanishing of the matched expression at the kinematical boundary 
y^ax) the further substitution [22] is made: 

L^L='-ln((^Y-( + A ' (2-17) 

P \\XLyJ \XLyma^J J 

where yo = 6 for y = C and yo = 1 otherwise. The values p = 1 and xl = 1 arc taken as 
default. The arbitrariness in the choice of the logarithm to be resummed can be quantified 
by varying the constant xl. 

3. Determination of the strong coupling constant 

As in the analysis of Ref. [27] we have studied the six event-shape distributions thrust 
T [17], heavy jet mass Mh [18], total and wide jet broadening {Bt,Bw) [19], C-parameter 
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C [20] and the two-to-three-jet transition parameter in the Durham algorithm — Inys [21]. 
The definitions of these variables and a discussion of their properties can be found in 
Refs. [1] and [22,50]. 

The measurements have been carried out by the ALEPH collaboration [1] ^ , at centre- 
of-mass energies of 91.2, 133, 161, 172, 183, 189, 200 and 206 GeV. Earlier measurements 
and complementary data sets from the LEP experiments and from SLD can be found in 
Refs. [2-6] . The event-shape distributions were computed using the reconstructed momenta 
and energies of charged and neutral particles. The measurements have been corrected for 
detector effects, i.e. the final distributions correspond to the so-called particle (or hadron) 
level. The particle level is defined by stable hadrons with a Ufetime longer than 10~^ s after 
hadronisation and Icptons according to the definition given in [51]. In addition, at LEP2 
energies above the Z peak they were corrected for initial-state radiation effects. At energies 
above 133 GeV, backgrounds from 4-fermion processes, mainly from W-pair production and 
also ZZ and Z7*, were subtracted following the procedure given in [1]. The experimental 
uncertainties were estimated by varying event and particle selection cuts. They are below 
1% at LEPl and slightly larger at LEP2. For further details we refer to Ref. [1]. 

The determination of the coupling constant from these data follows very closely the 
approach chosen in Refs. [1,27]. The perturbative predictions for the distributions, as 
described in section 2, arc calculated to the same order of perturbation theory for all of 
these variables and fit to the data. The measurements from several variables are combined, 
since this yields a better estimator of than using a single variable. Furthermore, the 
spread of values of is an independent estimate of the theoretical uncertainty. At centre- 
of-mass energies above the Z peak the statistical uncertainties are larger and background 
conditions are more difficult than at the peak. Therefore a combination of measurements 
is particularly important for those energies. We apply the same combination procedure as 
described in [1,27], which is based on weighted averages and takes into account correlations 
between the event-shape variables. However, in this paper we also investigate a combination 
procedure which excludes the perturbative uncertainty as weight, in order to evaluate the 
stability of our nominal combination method. Note that for energies above Mz we adopt 
the same treatment of statistical uncertainties as in [1,27]. The same holds for the fit 
ranges chosen at the LEP2 energies [27], whereas for the LEPl data the fit ranges have 
been slightly extended, motivated by the expected better description of the two-jet range 
by the resummed calculations. 

In this paper we present fits of matched NNLO-I-NLLA predictions and compare them 
to pure NNLO and matched NLO-I-NLLA calculations as used in the analyses of Refs. [1,27] . 
The nominal value for the renormalisation scale = fi/Q is unity. The perturbative QCD 
prediction is corrected for hadronisation and resonance decays by means of a transition 
matrix, which is computed with the Monte Carlo generators PYTHIA [52], HERWIG [53] 
and ARIADNE [54], all tuned to global hadronic observables at Mz [51]. The parton level 
is defined by the quarks and gluons present at the end of the parton shower in PYTHIA 
and HERWIG and the partons resulting from the colour dipole radiation in ARIADNE. 

^The tables with numbers and uncertainties for all variables can be found at 
http : / /aleph . web . cern . ch/aleph/ QCD/alephqcd . html. 
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Corrected measurements of event-shape distributions are compared to the theoretical calcu- 
lation at particle level. In section 6 we investigate the use of the NLO+LL event generator 
HERWIG++. 

The value of is determined at each energy using a binned least-squares fit. The fit 
programs of Ref. [27] have been extended to incorporate the NNLO-I-NLLA calculations. 
Only statistical uncertainties arising from the limited number of observed events, from the 
number of simulated events used to calculate hadronisation and detector corrections and 
from the integration procedure of the NNLO coefficient functions are included in the of 
the fit. Its quality is good for all variables at all energies. Nominal results for thrust and 
— lnj/3, based on (2.9) and using the fitted values of a^, are shown in Fig. 1, together with 
the measured distributions. The resulting measurements of ag{Q) for all six event shapes 
are given in Table 1 for 91.2 to 172 GeV and in Table 2 for 183 to 206 GeV. Comparisons 
of fits using different perturbative approximations are shown in Fig. 2 for the variables 
thrust and ys and results for all variables at LEPl are given in Table 3. 




0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 t 1 « 2 4 6 8 10 



Figure 1: Distributions measured by ALEPH, after correction for backgrounds and detector 
effects, of thrust and the two-to-three-jet transition parameter in the Durham algorithm at energies 
between 91.2 and 206 GeV, together with the fitted NNLO-^NLLA QCD predictions. The error 
bars correspond to the statistical uncertainties. The fit ranges cover the central regions indicated 
by the soUd curves, the theoretical predictions extrapolate well outside these fit ranges, as shown 
by the dotted curves. The plotted distributions are scaled by arbitrary factors for presentation. 

4. Systematic Uncertainties of ocg 

For a description of the determination and treatment of experimental systematic uncer- 
tainty we refer to Refs. [1, 27], since the identical approach is taken for this analysis. 
Similarly, the analysis of theoretical uncertainties goes along the lines of these earlier pub- 
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Figure 2: Distributions measured by ALEPH at LEPl, after correction for detector effects, of 
thrust and the two-to-three-jet transition parameter in the Durham algorithm. Fitted QCD pre- 
dictions at different orders of perturbation theory are overlaid. The lower insets show a relative 
comparison of data and QCD fits. 

lications. The main source of arbitrariness in the predictions is the choice of the renormal- 
isation scale and of the logarithmic rescaling variable xl- The residual dependence of 
the fitted value of as{Mz) on the renormalisation scale is shown in Fig. 3, for the same 
two variables as in the previous figures. Most notably, the matching of NLLA terms to the 
NNLO prediction does not lead to a reduced scale dependence, compared to pure NNLO 
only, but at least to an improvement compared to NLO-I-NLLA. This could be anticipated 
by the discussion in section 2 on the scale dependence of the NNLO and NLLA predictions. 
A further study of this particular aspect is described in section 6 below. 

The systematic uncertainty related to missing higher orders is estimated with the 
uncertainty-band method recommended in Ref. [22]. Briefly, this method derives the 
uncertainty of from the uncertainty of the theoretical prediction for the event-shape 
distribution and proceeds in three steps. First a reference perturbative prediction, here 
NNLO-I-NLLA with = 1 and = 1, is determined using the value of ag obtained from 
the combination of the six variables and eight energies, as explained in section 5. Then 
variants of the prediction with different choices for and x l , for the kinematic constraint 
ymax and the modification degree power p are calculated with the same value of a^. A 
variation of the matching scheme as advocated in Ref. [22] was not included in the list 
of variants, since no i?-matching scheme is presently available at NNLO+NLLA. In each 
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Figure 3: Dependence of the extracted a^-value on the renormahsation scale when fitting the 
distributions of thrust (left) and the two-to-three-jet transition parameter in the Durham algorithm 
(right) with predictions at NNLO+NLLA (sohd), NNLO (dashed) and NLO+NLLA (dotted). 

bin of the distribution for a given variable, the largest upward and downward differences 
with respect to the reference prediction are taken to define an uncertainty band around 
the reference theory. In the last step, the value of in the reference prediction is varied, 
in order to find the range of values which result in predictions lying inside the uncertainty 
band for the fit range under consideration. In contrast to the original method [22] we do 
not require the reference prediction to lie strictly inside the uncertainty band, since for 
the present NNLO-I-NLLA calculations the latter is still subject to statistical fluctuations. 
Instead, we make a fit of the reference theory with as free parameter to the uncertainty 
band, which includes the statistical uncertainty on the C coefficient, as in Ref. [27]. The 
values of fitted to the upper and lower contour of the uncertainty band finally set the 
perturbative systematic uncertainty. The upward and downward uncertainties are very 
similar in magnitude and the larger is quoted as symmetric uncertainty. The method is 
illustrated in Fig. 4 for thrust and the two-to-three-jet transition parameter in the Durham 
algorithm —Iny^. 

The combined value of a^, used to derive the systematic perturbative error, depends 
itself on the theoretical error. Hence the procedure of calculating the combination and 
its perturbative error is iterated until convergence is reached, typically after two iterations. 

At LEP2 energies the statistical fluctuations are large. In order to avoid biases from 
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Figure 4: Theoretical uncertainties for the distributions of thrust (left) and the two-to-three-jet 
transition parameter in the Durham algorithm (right) at LEPl. The filled area represents the 
perturbative uncertainties of the distribution for a given value a^. The curves show the reference 
prediction with ± Aa <, . The theoretical uncertainty Aa<, is derived from a fit of the reference 
theory to the contour of the uncertainty band for the actual fit range. 

downward fluctuations, the theoretical uncertainties are calculated with the value of 
obtained by the global combination procedure. For each energy point and in each variable, 
the combined is evolved to the appropriate energy scale and the uncertainty is calculated 
for the fit range used for the different variables. 

An additional error is evaluated for the b-quark mass correction procedure. This 
correction has only been calculated to 0(0^) for the differential coefficients; no resummed 
and NNLO expressions are yet available. We have updated the calculations for the massive 
coefficients used in [1,27] to include now three different sets for Mb = 4.0, 4.5 and 5.0 
GeV/c^ . The difference in obtained with these different sets is taken as systematic 
error. The difference between the massless and massive expression for the hadronic cross 
section is already rather small and not included in this estimate. 

The total perturbative uncertainty quoted in the tables is the quadratic sum of the 
errors for missing higher orders and for the mass correction procedure. The total pertur- 
bative error is between 3% and 5% at Mz and decreases to between 2% and 3% at LEP2 
energies. 

The hadronisation model uncertainty is estimated by comparing the standard hadron- 
level event generator programs HERWIG and ARIADNE to PYTHIA for both hadroni- 



- 12 - 



sation and detector corrections. The same set of corrections as in Ref. [1] is used. Both 
corrections are calculated with the same generator in order to obtain a coherent descrip- 
tion at the hadron level. The maximum change with respect to the nominal result using 
PYTHIA is taken as the systematic error. At LEP2 energies the hadronisation model 
uncertainty is again subject to statistical fluctuations. These fluctuations are observed 
from one energy to the next and originate from limited statistics of the fully simulated 
detector-correction functions. Since non-perturbative effects are expected to decrease with 
1/(5, the energy evolution of hadronisation errors has been fitted to a simple A + B/Q 
parametrisation. The fit was performed for each variable separately. In the fit procedure a 
weight scaling with luminosity is assigned to the hadronisation uncertainty at each energy 
point. This ensures that the hadronisation uncertainty at Mi, which is basically free of 
statistical fluctuations, is not altered by the procedure. As in the case of experimental 
systematic uncertainties, the hadronisation uncertainty is essentially identical to that pub- 
lished in [1] . In section 6 we present an attempt to use modern event generators to estimate 
the hadronisation corrections. 

The perturbative component of the error, which is the dominant source of uncertainty 
in most cases, is highly correlated between the energy points. The perturbative errors 
decrease with increasing Q, and faster than the coupling constant itself. The overall error 
is in general dominated by the renormalisation scale dependence. 

5. Combined Results 

The measurements obtained from the six different variables using NNLO-I-NLLA calcula- 
tions are combined into a single measurement per energy using weighted averages. The 

same procedure as in Refs. [1,27] is applied here. However, we investigate also the impact 
of the theoretical uncertainties in the calculations of the weights, as described in section 6 
below. 

In Table 7 the weighted averages are given for all LEPl and LEP2 data, as well as for 
the LEP2 data only. Essentially identical results with very similar errors are found. The 
fitted values of at the various centre-of mass energies are displayed in Fig. 5 and com- 
pared to the QCD three-loop formula for the running of the coupling constant. Excellent 
agreement of the data with the expected energy dependence is observed. 

In Table 8 we compare the combined results obtained at NNLO+NLLA accuracy to 
the results at NNLO and NLO-I-NLLA. The numbers given in Table 8 supersede those 
published in [27] at NNLO and NLO-I-NLLA and can be traced back to the following 
changes in the present analysis 

• for the normalisation to cThad an expansion of the ratio c7o/<Thad was applied in Refs. [1, 
27], while here the exact value is used; 

• new massive coefficients using Mb = 4.5 GeV /(? up to NLO are used; 

• a massive expression for Uhad is iiow adopted; 
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Figure 5: The measurements of the strong couphng constant as between 91.2 and 206 GeV. 
The results using the six different event-shape variables are combined with correlations taken into 
account. The inner error bars exclude the perturbative uncertainty, which is expected to be highly 
correlated between the measurements. The outer error bars indicate the total error. A fit of the 
three-loop evolution formula using the uncorrelated errors is shown. The shaded area corresponds 
to the uncertainty in the fit parameter = 284 ± 14 MeV of the three-loop formula, eq. (2.3). 



• for a given observable and energy, the same fit ranges (given in Tables 1 and 2) are 
applied to different theoretical predictions; 

• a small transcription error in the fit program used in [27] when calculating the NNLO 
term for —Iny^ is corrected; 

• the previously incomplete treatment of large-angle soft radiation [14] in EERAD3 is 
corrected, resulting in minor numerical shifts in the NNLO coefficients; 

• while in [1,27] the NL04-NLLA predictions were obtained by a numerical derivative 
of R (cf. eq. (2.5)), we now compute the differential distributions analytically for the 
resummed part, which yields a better numerical stability. We apply this procedure 
also to NNLO-FNLLA. 

As was anticipated in section 2, the matching of the NNLO prediction with the re- 
summation at NLLA introduces a renormalisation scale dependence which is absent in the 
pure NNLO C3(S6, clS described in detail in section 2. This is reflected by the increased 
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perturbative and finally total uncertainty of the NNLO+NLLA result compared to NNLO, 
as can be seen by comparing Table 5 for the combined value of as{Mz) at different energies 
at NNLO+NLLA with Table 6 at NNLO. However, compared to the NLO+NLLA fit, an 
improvement of more than 20% is obtained for the perturbative error. The central values 
of the fits for the different approximations turn out to be pretty similar. The fitted values 
of the coupling constant as found from the various event-shape variables, combined over 
all energies, are shown in Fig. 6. Besides the larger uncertainties, at NNLO+NLLA we 
observe the same reduced scatter of the results compared to NLO+NLLA as already re- 
ported previously [27]. However, the effect is not as strong as going from a NLO fit (where 
the scatter is largest) to a pure NNLO fit. 
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Figure 6: The measurements of tlie strong coupling constant as for the six event shapes, at 
y/s — Mz, when using QCD predictions at different approximations in perturbation theory. The 
shaded area corresponds to the total uncertainty, as in Fig. 5. 



6. Systematic studies 

6.1 lni?(;u)-matching scheme 

As described in section 2, we have computed the two-loop terms proportional to the renor- 
malisation scale in the resummation and matching functions (eq. 2.14) and recomputed 
the theoretical error in the new matching scheme, which we call the lni2(/x)-scheme. It 
is important to note that this new matching scheme does not affect the central values of 
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the individual fit results, since the scale = xl = 1 is used. Only the perturbative 
uncertainty will be changed because of the different scale dependence. However, since this 
uncertainty enters as a weight in the combination procedure, and different event shapes 
display different scale dependence with the lni2(/i)-matching scheme, the central values 
also change in the combined results. The results for the LEPl centre-of-mass energy are 
given in Table 9 for all six variables, whereas the combination of all variables and energies 
is listed in Table 10. The corresponding uncertainty bands are shown in Fig. 7 for thrust 
and the three-jet transition variable. It can be seen that in this modified matching scheme 
the renormalisation scale and xl dependence are very strongly reduced, leading to a more 
precise determination. However, given the fact that for a consistent analysis the full 
NNLLA calculation should be matched to the NNLO prediction, we prefer to quote the 
values obtained with the standard lni?-matching as our main result. 




Figure 7: Theoretical uncertainties for the distributions of thrust (left) and the three-jet transition 
variable (right) at LEPl, using the lni?(/i)-niatching scheme at NNLO+NLLA. 

6.2 Normalisation and quark mass effects 

In our nominal analysis the theoretical prediction is normalised to the total hadronic cross 
section, taking properly into account the production of massive b-quarks. Furthermore, 
mass corrections are applied for the fixed-order coefficients at leading and next-to-leading 
order. In order to study the impact of different normalisation and mass correction schemes 
the analysis has been repeated with alternative options, as summarised in Table 11 for the 
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LEPl data. The observed differences when using either the massive or massless hadronic 
cross section as normahsation are rather minor (first and second row in Table 11). The 
alternative approach to applying the exact correction cxo / crhad j namely expanding this ratio 
and correspondingly changing the fixed order coefficient functions B ^ B and C ^ C, 
has been adopted in previous publications. In this case the results are lowered by about 
0.5% (third row in Table 11). For completeness, we also give the results obtained with 
massless coefficients throughout (fourth row in Table 11). This again lowers the result for 
as(Mz) by 0.5% at LEPl, but has almost no effect at LEP2. The last two rows give the 
results obtained with different values for the b-quark pole mass, which we use to derive the 
uncertainty for the mass correction procedure. 

6.3 Combination method 

Our nominal combination procedure is based on weighted averages with weights propor- 
tional to the total significance i.e. oc l/a^g^, where utot is the total uncertainty for an 
individual measurement, thus including the perturbative error in the weight calculation. 
However, it can be expected that the theoretical uncertainty for a given observable is highly 
correlated between different energies, the main de-correlation effect being related to the 
different fit ranges. In contrast, the theoretical uncertainties of different variables at the 
same energy are clearly less correlated, since missing higher-order contributions are likely 
to be different. Therefore it is instructive to study the stability of the combination method 
by using only the largely uncorrelated uncertainty component as weight when combining 
the results from different energies, while for the first-step combination of different variables 
at the same energy the theoretical uncertainties are still included. As a result of such a 
procedure the importance of the statistical error is significantly enhanced, leading to a 
reduced weight of the LEP2 data with respect to LEPl. 

In Table 12, the newly obtained weights are compared to the nominal weights for 
the combination of measurements at different energies, and the resulting combined mea- 
surements are given in Table 13. As anticipated, now the overall combination as well as 
the resulting theoretical uncertainty are dominated by the LEPl data, while almost no 
difference is observed when combining only the LEP2 energy points. 

6.4 Hadronisation corrections from NLO-|-LL event generators 

In recent years substantial progress has been achieved in the development of modern Monte 
Carlo event generators targeted in particular towards the LHC era and often implemented 
in object oriented C++ frameworks. Compared to the legacy generators used in the LEP 
era, the new programs include in part NLO corrections matched to parton showers at 
leading logarithmic accuracy (LL) for various processes. Here we use HERWIG++ [33] 
version 2.3 for our investigations, which is based on ThePEG [55], a general framework for 
implementing Monte Carlo generator classes. The nominal version for HERWIG++ uses 
a LO+LL configuration which features matrix element corrections for the matching of the 
hard scattering process to the parton showers. Furthermore, two schemes for the imple- 
mentation of NLO corrections, namely the MCNLO [56] and POWHEG [57] schemes, are 
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available^. The actual implementations of the general NLO to LL matching prescriptions 
are given in Ref. [58] for MCNLO and in Ref. [59] for POWHEG. Technically, a simula- 
tion at NLO+LL is obtained in two steps, where first the NLO partonic configurations are 
generated [60] and second these events are passed to HERWIG++ to simulate the parton 
shower, hadronisation and resonance decays. 

It should be noted that the nominal LO+LL version with matrix element corrections 
of HERWIG++ has been extensively tuned to a variety of experimental data, including jet 
rates, event shapes and particle multiplicities in e+e^ annihilations from LEP and heavy 
flavour data from the B-factories [61], in order to obtain the best possible set of parameters. 
However, according to [61], the quality of this fit, with an overall per degree of freedom 
between 5 and 6, is limited, and this is, to some extent, related to a slightly overestimated 
amount of gliion radiation in the parton shower. It can not be expected that this tuning 
is optimal for the NLO+LL versions of HERWIG++. A re-tuning of the parameters 
using MCNLO and POWHEG is beyond the scope of this paper, but we used parameters 
suggested in Ref. [62] for MCNLO and checked a few of the main parameters for POWHEG 
using the ALEPH event shapes [1]. To this end four parameters were individually varied 
in a simplified grid-search procedure in which for each configuration the with respect 
to the full range of the six event-shape distributions from ALEPH at LEPl was recorded. 
This resulted in an improved description of the distributions studied here, presumably at 
the expense of other distributions included in the more general tuning of HERWIG-I— f. 

The following parameters were determined according to this procedure for POWHEG: 
AlphaMZ = 0.134, cutoftKinScale = 2.7, PSplitLight = 1.1 and PwtDIquark = 
0.6. The meaning of these parameters can be inferred from [33]. In addition, the partonic 
configuration was generated with a value of A = 170 MeV to set the scale for the hardest 
gluon emission according to the evolution in the MS-schcmc, yielding as(Mz)=0.11. 

We compare in a first step the prediction for the event shape distributions of HER- 
WIG++ to both the high precision data at LEPl from ALEPH and the predictions from 
the legacy generators PYTHIA, HERWIG and ARIADNE. We recall that the latter have 
all been tuned to the same global QCD observables measured by ALEPH [51] at LEPl, 
which included cvcnt-sliapc variables similar to the ones analysed here. In Fig. 8 the gen- 
erator predictions for thrust, — ln(2/3) and the total and wide jet broadenings are compared 
to the ALEPH data. 

In general it appears that the shape of HERWIG++ is similar to both HERWIG-f-f- 
with MCNLO and HERWIG, but all differ in normalisation. A better description is ob- 
tained using HERWIG-F-F with POWHEG. PYTHIA and ARIADNE yield by far the best 
description. To quantify the performance of the generators, in Table 14 we have compiled 
the x^ of their predictions with respect to the event shapes studied at LEPl, including the 
experimental systematic uncertainty. A complete re-tuning of the HERWIG++ parameters 
to the same data used to tune PYTHIA, HERWIG and ARIADNE would likely improve 
their performance. 

To investigate the origin of the observed differences between the generators, it is in- 

^We use the notation MCNLO for the method, while MC@NLO denotes the program. 
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Figure 8: Residuals of liadron level Monte Carlo predictions with respect to the ALEPH data. 
The shaded area indicates the experimental uncertainty. 



structive to consider the parton-level predictions and the hadronisation corrections sepa- 
rately. The parton level predictions from the generators are calculated with final state par- 
tons at the end of the parton shower. These are compared to the complete NNLO-I-NLLA 
calculation in Fig. 9. For thrust and in particular the total jet broadening a reasonable 
agreement between NNLO+NLLA and HERWIG-h+ with POWHEG, as well as a fair 
agreement with PYTHIA and ARIADNE is observed, while other HERWIG variants show 
a clear deviation. For — ln(y3) and the wide jet broadening all legacy generators provide a 
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satisfactory description and HERWIG++ based predictions exhibit some systematic differ- 
ences in sliape. It is worth noting that for thrust and the total jet broadening, the PYTHIA 
prediction overestimates the NNLO+NLLA calculation by about 10%, with the shape in 
reasonable agreement, whereas a better agreement is seen for the wide jet broadening and 
-ln(y3)- 
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Figure 9: Normalised residuals with respect to the NNLO+NLLA calculation of parton level 
predictions obtained from different Monte Carlo generators. 

The hadronisation corrections to be used in the fits to the data arc shown in Fig. 10. 
HERWIG++ with POWHEG yields a similar shape as the legacy programs, but differs 
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in the normalisation. The other HERWIG++ predictions differ most notably in shape 
from the former. In Table 15 the fit results obtained with all generators for hadronisation 
corrections are given. In most cases, the fits based on HERWIG++ and HERWIG++ 
with MCNLO are significantly worse than for the other generators, but for individual 
variables like the wide jet broadening an opposite behaviour is observed. The fit quality of 
HERWIG++ with POWHEG is similar to the outcome of the legacy generators. Given the 
similar shape but different normalisation of HERWIG++ with POWHEG, the resulting 
values of are significantly lower, overall by 3 %. 

7. Discussion and Conclusions 

We have performed a determination of the strong coupling constant from event-shape 
data measured by the ALEPH collaboration [1], based on the perturbative QCD results 
at next-to-next-to-leading order (NNLO) matched to resummation in the next-to-leading- 
logarithmic approximation (NLLA) [24]. 

Comparing our results to both the fit using purely fixed-order NNLO predictions [27] 
and the fits based on earlier NLLA-I-NLO calculations [1], we make the following observa- 
tions: 

• The central value obtained by combining the results for six event-shape variables and 
the LEPl and LEP2 centre-of-mass energies, 

as(Mz) = 0.1224 ± 0.0009 (stat) ± 0.0009 (exp) ± 0.0012 (had) ± 0.0035 (theo), 

is slightly lower than the central value of 0.1228 obtained from fixed-order NNLO 
only, and slightly larger than the NLO-I-NLLA results. We note that in this analysis 
an improved normalisation to the total hadronic cross section has been used, which 
leads to minor deviations to previously reported results. 

The fact that the central value is almost identical to the purely fixed-order NNLO 
result could be anticipated from the findings in Rcf. [24]. There it is shown that in the 
three-jet region, which provides the bulk of the fit range, the matched NLLA-I-NNLO 
prediction is very close to the fixed-order NNLO calculation. 

• The dominant theoretical uncertainty on as{Mz), as estimated from scale variations, 
is reduced by 20% compared to NLO-hNLLA. However, compared to the fit based 
on purely fixed-order NNLO predictions, the perturbative uncertainty is increased in 
the NNLO-I-NLLA fit. The reason is that in the two-jet region the NLLA-I-NLO and 
NLLA-I-NNLO predictions agree by construction, because the matching suppresses 
any fixed order terms. Therefore, the renormalisation scale uncertainty is dominated 
by the ncxt-to-lcading-logarithmic approximation in this region, which results in a 
larger overall scale uncertainty in the fit. 

• As already observed for the fixed-order NNLO results, the scatter among the values 
of as{Mz) extracted from the six different event-shape variables is smaller than in 
the NLO+NLLA case. 
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Figure 10: Hadronisation corrections (ratio of hadron to parton level predictions) obtained from 
different Monte Carlo generators. 



• The matching of NLLA+NNLO introduces a mismatch in the cancellation of renor- 
maUsation scale logarithms, since the NNLO expansion fully compensates the renor- 
malisation scale dependence up to two loops, while NLLA only compensates it up to 
one loop. In order to assess the impact of this mismatch, we introduced the lni?(/i) 

matching scheme, which retains the two-loop renormalisation terms in the resummed 
expressions and the matching coefficients. In this scheme, a substantial reduction of 
the perturbative uncertainty from ±0.0035 (obtained in the default lni?-scheme) to 
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±0.0022 is observed, which might indicate the size of the ultimately reachable pre- 
cision for a complete NNLO+NNLLA calculation including the currently unknown 
resummed function 33 for all shape variables. Although both schemes are in principle 
on the same theoretical footing, it is the more conservative error estimate obtained 
in the lni2-scheme which should be taken as the nominal value, since it measures the 
potential impact of the yet uncalculated finite NNLLA-terms. 

• Bottom quark mass effects, which are numerically significant mainly at the LEPl 
energy, were included through to NLO. Compared to a purely massless evaluation 
of the distributions, the inclusion of these mass effects enhances as{Mz) by 0.8%. 
Compared to the previously used expansion of the mass corrections, an enhancement 
of 0.4% is observed. 

• The averaging of as{Mz) values obtained at the various LEPl and LEP2 energies 
weights the different measurements by their total uncertainties. Excluding the error 
on the perturbative prediction from this weighting enhances the importance of the 
very precise LEPl data over the LEP2 data, and yields an as{Mz) value which is 
lower than our default result by only 0.2%, thereby demonstrating the very good 
consistency of the LEPl and LEP2 results. 

• We have investigated hadronisation corrections obtained from NLO+LL parton shower 
simulation using HERWIG++ with two different schemes. Results for based on 
corrections from HERWIG++ with POWHEG are slightly lower than with nomi- 
nal corrections from PYTHIA. Comparing hadron level predictions with data reveals 
that HERWIG++ with POWHEG yields an improved prediction over HERWIG, 
HERWIG++ and HERWIG++ with MCNLO, but docs not reach the same level of 
agreement as PYTHIA and ARIADNE. Further, we observe a certain discrepancy 
between MCNLO and POWHEG, which might indicate unresolved tuning issues. 
Therefore, while the first studies with HERWIG++ look rather promising, we retain 
for the time being PYTHIA as generator for our nominal result. 

• From the study of hadronisation corrections we also make the following important 
observation. It appears that there are two "classes" of variables. The first class 
contains thrust, C-parameter and total jet broadening, whereas the second class 
consists of the heavy jet mass, wide jet broadening and the two-to-three-jet transition 
parameter — Inya. For the first class, using the standard hadronisation corrections 
from PYTHIA, we obtain as{Mz) values around 0.125 — 0.127, some 5% higher than 
those found from the second class of variables. In a study of higher moments of event 
shapes [32], indications were found that variables from the first class still exhibit 
sizable missing higher order corrections, whereas the second class of observables have 
a better perturbative stability. In this paper, from Fig. 9, we observe that this 
first class of variables gives a parton level prediction with PYTHIA, which is about 
10% higher than the NNLO+NLLA prediction. The PYTHIA curve is obtained 
with tuned parameters, where the tuning to data had been performed at the hadron 
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level. Indeed, this tuning results in a rather large effective coupling in the parton 
shower, which might partly explain the larger parton level prediction of PYTHIA. 
However, since the tuning has been performed at hadron level, this implies that the 
hadronisation corrections come out to be smaller than what would have been found by 
tuning a hypothetical Monte Carlo prediction with a parton level corresponding to the 
NNLO+NLLA prediction. Thus, in the end, the PYTHIA hadronisation corrections, 
applied in the fit, might be too small, resulting in a larger as{Mz) value. Such 
problems do not appear to exist for the second class of variables. 

In summary, there are indications that the first class of event shapes still suffers 
from significant missing higher order contributions, even beyond NNLO+NLLA. This 
might also have led to a tuning of parton shower models which underestimates the 
hadronisation corrections for these variables, and consequently results in somewhat 
larger values of the fitted strong coupling. Since up to now the hadronisation un- 
certainties have been estimated from the differences of parton shower based models, 
tuned to the data, it is likely that for these event shapes the uncertainties were 
underestimated and not able to account for a possible systematic shift. 

In future work it would be interesting to investigate the effect of NNLLA resummation 
terms for all six event shapes, of electroweak corrections, of quark mass effects beyond NLO 
and of non-perturbative power-law corrections as well as further studies with HERWIG++, 
in particular using the newly developed improved algorithm for merging matrix elements 
with angular-ordered parton showers [63]. 
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Q = 91.2 GeV 


V (All XCAiPs-ZAv^ 


T 






c 








0.1265 


0.1186 


0.1211 


0.1252 


0.1196 


0.1268 


stat. error 


0.0002 


0.0002 


0.0003 


0.0002 


0.0002 


0.0002 


exp. error 


0.0008 


0.0011 


0.0010 


0.0007 


0.0007 


0.0007 


pert, error 


0.0048 


0.0029 


0.0033 


0.0050 


0.0046 


0.0053 


hadr. error 


0.0019 


0.0017 


0.0042 


0.0016 


0.0017 


0.0022 


total error 


0.0052 


0.0036 


0.0055 


0.0053 


0.0049 


0.0058 


fit range 


0.75-0.91 


1.6-4.0 


0.10-0.22 


0.36-0.74 


0.09-0.19 


0.16-0.30 


Q = 133 GeV 


variable 


T 




Mh 


C 


Bw 


Bt 


as 


0.1193 


0.1199 


0.1149 


0.1174 


0.1158 


0.1205 


stat. error 


0.0047 


0.0057 


0.0053 


0.0037 


0.0032 


0.0031 


exp. error 


0.0006 


0.0004 


0.0012 


0.0008 


0.0008 


0.0010 


pert, error 


0.0039 


0.0024 


0.0027 


0.0040 


0.0038 


0.0044 


hadr. error 


0.0015 


0.0010 


0.0027 


0.0012 


0.0010 


0.0013 


total error 


0.0063 


0.0063 


0.0067 


0.0057 


0.0051 


0.0056 


fit range 


0.75-0.94 


1.6-4.4 


0.08-0.25 


0.30-0.75 


0.08-0.25 


0.13-0.35 


Q = 161 GcV 


variable 


T 


-Inys 


Mh 


C 


Bw 


Bt 


Ois 


0.1172 


0.1183 


0.1225 


0.1190 


0.1186 


0.1238 


stat. error 


0.0080 


0.0082 


0.0072 


0.0066 


0.0047 


0.0052 


exp. error 


0.0006 


0.0004 


0.0012 


0.0008 


0.0008 


0.0010 


pert, error 


0.0036 


0.0022 


0.0025 


0.0037 


0.0035 


0.0040 


hadr. error 


0.0014 


0.0007 


0.0022 


0.0011 


0.0008 


0.0010 


total error 


0.0088 


0.0085 


0.0079 


0.0076 


0.0060 


0.0067 


fit range 


0.75-0.94 


1.6-4.4 


0.08-0.25 


0.30-0.75 


0.08-0.25 


0.13-0.35 


Q = 172 GeV 


variable 


T 




Mh 


C 


Bw 


Bx 


as 


0.1120 


0.1095 


0.1079 


0.1093 


0.1036 


0.1108 


stat. error 


0.0077 


0.0098 


0.0085 


0.0063 


0.0063 


0.0069 


exp. error 


0.0006 


0.0006 


0.0012 


0.0008 


0.0008 


0.0012 


pert, error 


0.0035 


0.0021 


0.0024 


0.0035 


0.0033 


0.0039 


hadr. error 


0.0013 


0.0006 


0.0020 


0.0010 


0.0007 


0.0010 


total error 


0.0085 


0.0100 


0.0091 


0.0074 


0.0072 


0.0081 


fit range 


0.75-0.94 


1.6-4.4 


0.08-0.25 


0.22-0.75 


0.08-0.25 


0.11-0.35 



Table 1: Results for as{Q) as obtained from NNLO-I-NLLA fits to distributions of event-shape 
variables at Q = = 91.2, 133, 161 and 172 GeV. 
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Q = 183 GeV 


variable 


T 


— InwQ 




c 




1 




0.1131 


0.1083 


0.1129 


0.1094 


0.1091 


0.1148 


stat. error 


0.0036 


0.0050 


0.0038 


0.0032 


0.0027 


0.0030 


cxp. error 


0.0007 


0.0007 


0.0012 


0.0011 


0.0008 


0.0011 


pert, error 


0.0034 


0.0021 


0.0023 


0.0034 


0.0033 


0.0037 


hadr. error 


0.0013 


0.0005 


0.0018 


0.0010 


0.0007 


0.0010 


total error 


0.0051 


0.0055 


0.0050 


0.0049 


0.0043 


0.0050 


fit range 


0.80-0.96 


2.4-4.8 


0.06-0.20 


0.22-0.60 


0.065-0.20 


0.11-0.30 


Q = 189 GeV 


variable 


T 




Mh 


C 


Bw 


Bx 


as 


0.1119 


0.1087 


0.1087 


0.1121 


0.1056 


0.1137 


stat. error 


0.0026 


0.0031 


0.0032 


0.0018 


0.0019 


0.0019 


exp. error 


0.0007 


0.0005 


0.0017 


0.0009 


0.0009 


0.0012 


pert, error 


0.0033 


0.0020 


0.0023 


0.0034 


0.0032 


0.0037 


hadr. error 


0.0012 


0.0005 


0.0018 


0.0010 


0.0006 


0.0010 


total error 


0.0045 


0.0038 


0.0046 


0.0040 


0.0039 


0.0044 


fit range 


0.80-0.96 


2.4-4.8 


0.06-0.20 


0.22-0.60 


0.065-0.20 


0.11-0.30 


Q = 200 GeV 


variable 


T 




Mh 


C 


Bw 


Bx 


as 


0.1078 


0.1065 


0.1020 


0.1109 


0.1047 


0.1081 


stat. error 


0.0027 


0.0032 


0.0034 


0.0021 


0.0019 


0.0024 


exp. error 


0.0007 


0.0005 


0.0019 


0.0008 


0.0008 


0.0013 


pert, error 


0.0032 


0.0020 


0.0023 


0.0033 


0.0032 


0.0036 


hadr. error 


0.0012 


0.0005 


0.0016 


0.0009 


0.0006 


0.0010 


total error 


0.0045 


0.0038 


0.0048 


0.0041 


0.0039 


0.0046 


fit range 


0.80-0.96 


2.4-4.8 


0.06-0.20 


0.22-0.60 


0.065-0.20 


0.11-0.30 


Q = 206 GeV 


variable 


T 


-Ws 


Mh 


C 


Bw 


Bx 


as 


0.1084 


0.1040 


0.1076 


0.1076 


0.1051 


0.1089 


stat. error 


0.0025 


0.0032 


0.0024 


0.0019 


0.0017 


0.0020 


exp. error 


0.0007 


0.0005 


0.0012 


0.0008 


0.0008 


0.0011 


pert, error 


0.0032 


0.0020 


0.0022 


0.0033 


0.0031 


0.0035 


hadr. error 


0.0012 


0.0004 


0.0016 


0.0009 


0.0006 


0.0010 


total error 


0.0043 


0.0038 


0.0039 


0.0040 


0.0037 


0.0043 


fit range 


0.80-0.96 


2.4-4.8 


0.04-0.20 


0.22-0.60 


0.05-0.20 


0.11-0.30 



Table 2: Results for as{Q) as obtained from NNLO-I-NLLA fits to distributions of event-shape 
variables at Q = ^/s =183, 189, 200 and 206 GeV. 
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T 


C 


Mh 


Bw 


Bx 




NNLO+NLLA 


0.1266 


0.1252 


0.1211 


0.1196 


0.1268 


0.1186 




0.16 


0.47 


4.4 


4.4 


0.84 


1.89 


stat. error 


0.0002 


0.0002 


0.0003 


0.0002 


0.0002 


0.0002 


NLO+NLLA 


0.1282 


0.1244 


0.1180 


0.1161 


0.1290 


0.1187 




0.74 


1.88 


14.5 


19.6 


9.7 


4.7 


stat. error 


0.0002 


0.0002 


0.0003 


0.0002 


0.0002 


0.0002 


NNLO 


0.1275 


0.1273 


0.1248 


0.1242 


0.1279 


0.1192 




1.16 


1.08 


4.1 


2.74 


0.50 


1.17 


stat. error 


0.0002 


0.0002 


0.0004 


0.0002 


0.0002 


0.0003 


fit range 


0.75 - 0.91 


0.36 - 0.74 


0.10 - 0.22 


0.09 - 0.19 


0.16 - 0.30 


1.6-4.0 



Table 3: Fit results for as{Mz) using different predictions of perturbative QCD, with the renor- 
malisation scale fixed to = Mz- 



Q [GeV] 


91.2 


133 


161 


172 


183 


189 


200 


206 




0.1221 


0.1179 


0.1201 


0.1086 


0.1112 


0.1099 


0.1067 


0.1066 


stat. error 


0.0001 


0.0029 


0.0043 


0.0052 


0.0023 


0.0016 


0.0017 


0.0015 


exp. error 


0.0008 


0.0008 


0.0008 


0.0008 


0.0009 


0.0009 


0.0008 


0.0009 


pert, error 


0.0041 


0.0036 


0.0033 


0.0032 


0.0031 


0.0030 


0.0029 


0.0029 


hadr. error 


0.0018 


0.0012 


0.0010 


0.0010 


0.0009 


0.0008 


0.0008 


0.0008 


total error 


0.0045 


0.0049 


0.0056 


0.0062 


0.0040 


0.0036 


0.0036 


0.0034 


RMS 


0.0038 


0.0023 


0.0026 


0.0029 


0.0027 


0.0030 


0.0030 


0.0019 



Table 4: Combined results for as{Q) using NNLO+NLLA predictions. 



Q [GcV] 


91.2 


133 


161 


172 


183 


189 


200 


206 


as(Mz) 


0.1221 


0.1251 


0.1316 


0.1190 


0.1235 


0.1225 


0.1196 


0.1200 


stat. error 


0.0001 


0.0033 


0.0052 


0.0063 


0.0028 


0.0020 


0.0022 


0.0019 


exp. error 


0.0008 


0.0010 


0.0010 


0.0011 


0.0011 


0.0011 


0.0011 


0.0011 


pert, error 


0.0041 


0.0038 


0.0036 


0.0035 


0.0034 


0.0033 


0.0033 


0.0032 


hadr. error 


0.0018 


0.0014 


0.0012 


0.0011 


0.0011 


0.0010 


0.0010 


0.0010 


total error 


0.0045 


0.0053 


0.0065 


0.0074 


0.0047 


0.0042 


0.0042 


0.0041 


RMS 


0.0038 


0.0026 


0.0032 


0.0035 


0.0033 


0.0037 


0.0038 


0.0024 



Table 5: Combined results for as{Mz) using NNLO+NLLA predictions. 
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yi.z 


1 QQ 

loo 


1 1 

loi 


1 '70 


1 CQ 

loo 


1 on 
loy 


zUU 


zUD 




U.lzoy 


U.lz/U 


U.lolo 


n 1 1 no 


U.lzzD 


U.lzo4 


U.lzUU 


n 1 ono 
U.lzUz 


bXctx. error 




u.uuoo 


U.UUOl 


u.uuoo 


u.uuzo 


u.uuzu 


n 0091 
U.UUZl 


001 Q 
u.uuiy 


exp. error 


0.0009 


0.0009 


0.0009 


0.0010 


0.0010 


0.0010 


0.0010 


0.0010 


pert, error 


0.0030 


0.0030 


0.0028 


0.0028 


0.0027 


0.0026 


0.0025 


0.0025 


hadr. error 


0.0018 


0.0014 


0.0012 


0.0012 


0.0011 


0.0010 


0.0010 


0.0010 


total error 


0.0037 


0.0048 


0.0060 


0.0070 


0.0041 


0.0036 


0.0036 


0.0034 


RMS 


0.0036 


0.0014 


0.0043 


0.0019 


0.0027 


0.0027 


0.0034 


0.0024 



Table 6: Combined results for Qs(Mz) using NNLO predictions. 



data set 


LEPl + LEP2 


LEP2 


as(Mz) 


0.1224 


0.1224 


stat. error 


0.0009 


0.0011 


exp. error 


0.0009 


0.0010 


pert, error 


0.0035 


0.0034 


hadr. error 


0.0012 


0.0011 


total error 


0.0039 


0.0039 



Table 7: Weighted average of combined measurements for as{Mz) obtained at energies from 91.2 
GeV to 206 GeV and the average without the point at ^/s = Mz using NNLO+NLLA predictions. 



theory input 


NNLO+NLLA 


NNLO 


NLO+NLLA 




0.1224 


0.1228 


0.1215 


stat. error 


0.0009 


0.0008 


0.0010 


exp. error 


0.0009 


0.0009 


0.0009 


pert, error 


0.0035 


0.0027 


0.0053 


hadr. error 


0.0012 


0.0012 


0.0012 


total error 


0.0039 


0.0032 


0.0056 



Table 8: Comparison of combined results obtained with different theoretical predictions on as{Mz) 
using ALEPH data at energies from 91.2 GeV to 206 GeV. 



variable 


T 




Mh 


C 


Bw 


Bx 


InR 


0.0017 
0.0047 


0.0028 
0.0029 


0.0025 
0.0033 


0.0030 
0.0049 


0.0031 
0.0045 


0.0025 
0.0053 



Table 9: Comparison of the theoretical systematic uncertainties for the Ini? and lni?(/i) matching 
schemes. Only the uncertainty for missing higher orders as obtained from the uncertainty band 
method are included, using as(Mz)=0.1224. The total perturbative uncertainty also accounts for 
the mass corrections, the latter are the same for both matching schemes. 
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data set 


LEPl + LEP2 


LEP2 




0.1227 


0.1226 


stat. error 


0.0008 


0.0010 


exp. error 


0.0009 


0.0010 


pert, error 


0.0022 


0.0021 


hadr. error 


0.0012 


0.0011 


total error 


0.0028 


0.0028 



Table 10: Weighted average of combined measurements for q;s(Mz), obtained at energies from 
91.2 GeV to 206 GeV and without the point at ^/s = Mz, using in all cases the lni?(/i) matching 
scheme. 



variable 


T 




Mh 


C 


Bw 


Bx 


o"had (massive) 


0.1266 


0.1186 


0.1211 


0.1252 


0.1196 


0.1268 


fjhad (massless) 


0.1266 


0.1187 


0.1212 


0.1252 


0.1196 


0.1268 


massless expansion, massive A, B 


0.1260 


0.1183 


0.1208 


0.1247 


0.1192 


0.1262 


massless expansion, massless A, B 


0.1256 


0.1179 


0.1215 


0.1242 


0.1188 


0.1253 


Chad (massive, Mb = 4.0 GeV/c^ ) 


0.1264 


0.1185 


0.1212 


0.1251 


0.1195 


0.1267 


Chad (massive, Mb = 5.0 GeV/c^ ) 


0.1268 


0.1189 


0.1210 


0.1252 


0.1198 


0.1270 



Table 11: Results on as(Mz) from LEPl data using different normalisation and mass correction 
schemes. 



Q [GeV] 


91.2 133 161 172 183 189 200 206 


with pert. err. 
w/o pert. err. 


14.0 10.3 6.8 5.3 13.1 16.7 16.4 17.5 
80.0 2.5 2.6 2.5 2.8 3.1 3.2 3.3 



Table 12: Weights (in per cent) of the different centre-of-mass energy points in the global combi- 
nation, with and without the inclusion of theoretical uncertainties. 
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data set 


LEPl + LEP2 


LEP2 




0.1222 


0.1228 


stat. error 


0.0003 


0.0013 


exp. error 


0.0007 


0.0010 


pert, error 


0.0039 


0.0034 


hadr. error 


0.0017 


0.0011 


total error 


0.0044 


0.0040 



Table 13: Weighted average of the combined measurements for as{Mz), based on weights which 
do not include the theoretical uncertainty. 





T 


C 


Mh 


Bw 


Bj' 




global 


PYTHIA 


0.55 


1.05 


1.9 


2.5 


0.57 


1.31 


1.30 


ARIADNE 


0.44 


0.75 


0.60 


0.97 


0.58 


0.52 


0.65 


HERWIG 


6.9 


5.3 


9.4 


9.8 


4.4 


4.0 


6.6 


HW++ 


17.5 


16.5 


18.2 


13.4 


9.4 


5.6 


13.5 


HW++ MCNLO 


9.6 


15.9 


9.2 


10.5 


11.8 


5.5 


10.4 


POWIIEG 


3.9 


11.2 


8.5 


(3.9 


3.5 


2.3 


G.l 



Table 14: Comparison of hadron level predictions from various event generators to the ALEPH 
event-shape data. The Table shows the values normalised to the number of experimental bins, 
including statistical and experimental systematic uncertainties of the data. 



ocs{Mz) 


T 


C 


Mh 


Bw 


Bx 




PYTHIA 


0.1266 


0.1252 


0.1211 


0.1196 


0.1268 


0.1186 




0.16 


0.47 


4.4 


4.4 


0.84 


1.89 


ARIADNE 


0.1285 


0.1268 


0.1234 


0.1212 


0.1258 


0.1202 




0.96 


0.52 


2.5 


3.1 


2.15 


1.41 


HERWIG 


0.1256 


0.1242 


0.1253 


0.1203 


0.1258 


0.1203 




0.5 


0.65 


4.4 


2.0 


2.15 


0.8 


HW++ 


0.1242 


0.1228 


0.1299 


0.1212 


0.1238 


0.1168 


xVNdof 


6.6 


3.2 


3.3 


1.33 


2.65 


0.56 


HW++ MCNLO 


0.1234 


0.1220 


0.1292 


0.1220 


0.1232 


0.1175 


xVNdof 


10.7 


4.2 


2.2 


1.1 


5.7 


0.69 


HW++ POWHEG 


0.1189 


0.1179 


0.1236 


0.1169 


0.1224 


0.1142 


xVNdof 


1.46 


2.55 


3.8 


3.9 


1.54 


0.56 



Table 15: Fit results for as(Mz) using LEPl data and NLLO+NLLA but different hadronisation 
corrections. In all cases the same detector corrections, obtained from a full detector simulation 
using PYTHIA as generator is applied. The statistical errors are essentially unaltered compared to 
those in Table 3. 
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